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ABSTRACT 

A framework is introduced for coupling the evolution of galactic magnetic fields sustained 
by the mean-field dynamo with the formation and evolution of galaxies in cold dark matter 
cosmology. Estimates of the steady-state strength of the large-scale and turbulent magnetic 
fields from mean-field and fluctuation dynamo models are used together with galaxy proper¬ 
ties predicted by semi-analytic models of galaxy formation for a population of spiral galaxies. 
We find that the field strength is mostly controlled by the evolving gas content of the galaxies. 
Thus, because of the differences in the implementation of the star formation law, feedback 
from supernovae and ram-pressure stripping, each of the galaxy formation models considered 
predicts a distribution of field strengths with unique features. The most prominent of them is 
the difference in typical magnetic fields strengths obtained for the satellite and central galaxy 
populations as well as the typical strength of the large-scale magnetic field in galaxies of 
different mass. 
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1 INTRODUCTION 

Magnetic fields have been detected in all spiral galaxies when ob¬ 
served at sufficient sensitivity and resolution jBeck et al.||1996| 
Beck & Wielebinski 20131. The most informative observational 
tracers of magnetic fields are partially polarised synchrotron emis¬ 
sion and Faraday rotation. Polarised radio emission, together with 
significant Faraday rotation, indicate the presence of a component 
of the galactic magnetic field ordered on scales comparable to the 
resolution of the observation, of order a few hundred parsecs for 
nearby galaxies, and exceeding the integral turbulent scale which 
is of order 50-100pc. The random (turbulent) magnetic field b is 
usually stronger than the large-scale component B, so the degree of 
polarisation does not typically exceed 10-20 per cent. An average 


and from cosmic rays, as all of them are of comparable magnitude 


total field strength in nearby galaxies is B = (B~ -ffo 2 ) 1 / 2 : 
with b/B ~ 1-3 (Beck et al. 


9flG 


1996), ranging from about 4/rG in 
M31 ( |Fletcher et al.|2004) to 15 pG in M51 ^Fletcher et al,|2011| >. 

Magnetic fields contribute significantly to the structure and 
evolution of the interstellar medium (ISM) and the host galaxy. 
They affect the accretion of gas by dark matter haloes ( [Rodrigues] 
|et al.f TOlO). as well as the outflows and inflows in galaxies that 
have already formed (Birnboim 2009). Magnetic fields also influ¬ 
ence outflows since they affect the multi-phase structure of the in¬ 
terstellar gas as they confine hot gas bubbles produced by super¬ 
novae (|Ferriere et al,|199f]|Korpi et al.|1999||Shukurov et al,|2004| 
Hanayama & Tomisaka 2006). The magnetic contribution to the 

overall structure of galactic gas discs is at least as important as that 
from other sources of interstellar pressure, i.e. thermal, turbulent 


< ] Parker |[l~979l po emen 1987] |Boulares & Cox||1990||Fletcher &1 


|Shukurov|2001) . Half of the total interstellar pressure is thus due to 
non-thermal contributions, and, hence, magnetic fields directly af¬ 
fect the mean gas density. In turn, this affects significantly the star 
formation rate; it is perhaps surprising that the role of magnetic 
fields and cosmic rays in galaxy evolution has avoided attention 
for so long jBirnboim et al.||2015) . Magnetic fields also regulate 
star formation locally by controlling the collapse and fragmentation 
of molecular clouds ( |Mac Low|2009| |Peters et al.|2011[ [Crutcher] 
J2012). Magnetic fields contribute to interstellar gas dynamics not 
only directly but also by confining cosmic rays ( jBerezinskii et al.| 
(1990). The latter are effectively weightless and so are capable of 
driving galactic outflows (winds and fountains) [Breitschwerdt & 


Komossa|2000[ [Everett et al.|2008||Uhlig et al.|2012| Booth et al. 


2013) thus providing negative feedback on star formation in galac 
tic discs jVeilleux et al.|2005]|Putman et al.|2012) . 

Future radio telescopes such as the Square Kilometre ArrajQ 
will allow a dramatic increase in the sensitivity and angular reso¬ 
lution of radio observations to permit not only detailed studies of 
nearby galaxies but also reliable and extensive measurements of 
galactic magnetic fields at higher redshifts l |Gaensler et al.|2004) . 
As we demonstrate here, the evolution of galactic magnetic fields 
is sensitive to a variety of poorly constrained parameters in galactic 
evolution models and can thus provide a sensitive diagnostic of the 
physical processes involved in galaxy formation. 

The most successful theory for the production of magnetic 
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fields observed in galaxies is the turbulent dynamo (for reviews, 
see |Beck et al.|19 96; Shukurov 20071. Two types of dynamo are 
expected to be important: the fluctuation dynamo, which produces 
a small-scale, random magnetic field of coherence lengths com¬ 
parable to the turbulent scale in the plasma, and the mean-field 
dynamo, which produces a large-scale magnetic field ordered on 
scales greater than the turbulent scale. The details of the two dy¬ 
namos are briefly discussed in Section[3] 

There have been several attempts to understand the behaviour 
of large-scale magnetic fields in an evolving galaxy but they have 
been limited by the quality of galactic formation and evolution 
models available. |Beck et al.|[l994j > suggested that the seed mag¬ 
netic field for the large-scale dynamo can be produced in a pro¬ 
togalaxy or young galaxy by the fluctuation dynamo acting on a 
short time scale (see also Kulsrud et al. 1997, Arshakian et al. 2009 
|Schleicher et al.|2010HSur et al.|20 1 2[|Seifried et al.|20T4] >, and thus 
the observed galactic magnetic fields can be generated in 1-2 Gyr 
for typical galactic properties. The seed magnetic field for the fluc¬ 
tuation dynamo can be produced by many mechanisms ( |Durrer &| 
such as battery effects in the first generation of stars 
aL|2005| or in a rotating protogalaxy (|Mishustin &| 

w 

The only evolutionary effect included by |Beck et al.|[l994] l is 
a variation of the thickness of the galactic ionised gas layer with 
time. | Arshakian et al.| ( 2009 201 1} discuss dynamo action at vari¬ 
ous stages of galaxy formation at a qualitative level with reference 
to more recent models of galaxy formation and evolution. 

Here we quantitatively model the large-scale galactic mag¬ 
netic fields produced by a mean-field dynamo in the framework of 
specific hierarchical galaxy formation models, with the aim of pre¬ 
dicting the statistical properties of galactic magnetic fields in a large 
sample of galaxies at different redshifts. We employ semi-analytic 
models (SAMs) of galaxy formation based on simple physical pre¬ 
scriptions (for a review, see |Baugh|20Q6| . These models produce 
synthetic catalogues of galaxies that are able to reproduce a wide 
range of observables (e.g., the galaxy luminosity and stellar mass 
functions). Starting from the stellar and gaseous masses, disc sizes, 
star formation rates and circular velocities, obtained from SAMs 
for late-type galaxies, we calculate magnetic fields using dynamo 
theory. The galaxy formation models remain uncertain regarding 
the details of the physical processes involved and their relative im¬ 
portance. In particular, the nature and intensity of the feedback of 
star formation on the galactic discs, which is in turn sensitive to the 
form and strength of the interstellar magnetic field, remains unclear. 

We derive the dependence of the strength of magnetic field on 
the galactic mass, its evolution with redshift, and the statistical dis¬ 
tribution of magnetic field strengths in galactic samples. Remark¬ 
ably, different SAMs lead to distinct predictions regarding galac¬ 
tic magnetic fields. Thus, magnetic fields observed in high-redshift 
galaxies via their polarisation and Faraday rotation can help to re¬ 
fine galaxy formation models. 

An alternative approach to galaxy formation and evolution in a 
cosmological context involves the numerical solution of the magne¬ 
tohydrodynamic (MHD) equations for interstellar gas in the (evolv¬ 
ing or static) gravity field obtained from TV-body simulations. Gas 
dynamics is often simulated with adaptive mesh refinement (or a 
particle based technique such as Smooth Particle Hydrodynamics) 
in order to achieve higher spatial resolution in denser regions. In 
the best MHD simulations available now, the highest resolution is 
20-300pc ( |Wang & Abel|2009||Pakmor et al,|2014| . The driving of 
interstellar turbulence by supemovae is neglected or parameterised. 
Moreover, the turbulent dynamo action responsible for both large- 


Neronov]2013 


(Hanayama et 

Ruzmaikin 197 


scale and turbulent interstellar magnetic fields obviously needs a 
fully resolved turbulent flow to be modelled accurately. As a result, 
the processes responsible for the large-scale magnetic field cannot 
be captured at all, and the simpler fluctuation dynamo is controlled 
by random flows at unrealistic scales. Even so, such results should 
be treated with care as turbulent magnetic fields are known to be 
sensitive to flows well within the inertial range of interstellar turbu¬ 
lence at scales of a fraction of parsec ( [Brandenburg & Subramanianj 
2005). Thus, magnetohydrodynamic simulations of galaxy forma¬ 
tion and evolution need to be complemented with semi-analytic 
models coupled with a mean-field description of the large-scale 
magnetic field and models of the turbulent fields. The advantages 
of this approach include the opportunity to explore the parameter 
space and clarify the role of various physical effects. This oppor¬ 
tunity is particularly attractive given that many important physical 
processes are already heavily parameterised in the hydrodynamic 
and magnetohydrodynamic models. 

The paper is organised as follows. In Section [2] we discuss 
semi-analytic models of galaxy formation and evolution and Sec¬ 
tion [3] describes the generation of galactic magnetic fields by dy¬ 
namo mechanisms. Our results for the distribution and evolution of 
galactic magnetic fields can be found in Section [4] together with 
a discussion on the systematic connection between the gas content 
of galaxies and magnetic fields, the role of magnetic helicity dif¬ 
fusion, the redshift evolution of galactic magnetic fields emerging 
from our models and the dependence of magnetic field on the star 
formation rate. Our conclusions are presented in Section [5] The 
appendices provide further details of the models. 


2 GALACTIC EVOLUTION 


To obtain the galaxy properties needed to compute the strength of 
galactic magnetic fields we use semi-analytic models of galaxy for¬ 
mation (SAMs). These break down galaxy formation and evolution 
into a set of differential equations, each of which models a sepa¬ 
rate physical process, including: the merger history of dark matter 
haloes; the radiative cooling of gas inside these haloes and the sub¬ 
sequent formation of galactic discs; star formation in galactic discs; 
galaxy mergers driven by dynamical friction and bursts of star for¬ 
mation associated with them; the feedback due to supernovae and 
AGNs and the chemical enrichment of stars and gas. (For reviews 
on these topics, see |Baugh|2006||Benson|2010| ) 

We use three versions of the Durham SAM, GALFORM (Cole 
|etah|2000, Bower et al. 2006j, the version presented by|Baugh et al. 


( 2005 [ , hereafter referred to as BAU, that of|Lagos et al 1 2012 1 , re¬ 


ferred to as LAG and that of |Font et al.|[2008| , hereafter FON. All 
models reproduce the galactic luminosity function at redshift z = 0 
in the K and bj bands. The BAU version of GALFORM successfully 
reproduces the counts and redshift distribution of galaxies selected 
according to their sub-millimetre luminosity, and the abundances 
of Lyman-break galaxies [Lacey et al.|201l| and Lyman a emit¬ 
ters [Orsi et al. 2008). The FON version provides a better match to 
the observed colours of satellite galaxies. The LAG version repro¬ 
duces the atomic and molecular hydrogen mass functions l |Lagos| 
|et al.|20lT), the A-ba nd luminosity function at redshifts z = 1 and 
2 [Lagos et al.|2011| > and the CO(l-O) luminosity function. Both 
LAG and FON are derived from the model of B ower et al. (2006) 
and predict the properties of the AGN population [Fanidakis et al.| 

|20TT||20T2) . 

The choice of these models allows us to assess the variabil¬ 
ity of the predictions between state-of-art SAMs and the impact of 
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these differences on the the galactic magnetic fields inferred from BAU also assumes the Kennicutt IMF for quiescent star formation, 
them. but a top-heavy IMF is used for the starbursts. 


2.1 Semi-analytic models of galaxy formation and evolution 

The SAMs explored differ in various aspects, briefly described 
here. Further details can be found in the original papers referred 
to in the previous section. 


2.1.1 Cosmological model and dark matter merger trees 

The models differ in their assumed cosmologies and dark matter 
merger trees. The BAU model adopts a present-day cosmological 
constant energy density parameter £2 A =0.7, a matter density pa¬ 
rameter £2 m = 0.3, a Hubble parameter h = 0.7 and a spectrum 
normalisation Cg = 0.93. The LAG and FON models adopt the cos¬ 
mological parameters of the Millennium Simulation (Springel et al.| 
[2005) , £2 A = 0.75, £2 m = 0.25, h = 0.73 and <jg = 93. 

The dark matter (DM) merger trees used in the LAG and FON 
models are extracted from the Millennium IV-body cosmological 
simulation (Springel et al. 2005) , which has a halo mass resolution 
of 1.72 x 1O iu / 7 -1 M0. In the BAU model, the merger trees are built 
using the Monte-Carlo procedure described by |Cole et al.||2000) , 
with the minimum final halo mass for the merger trees at z = 0 set 
to 5 x 10 9 /j - 1 M q . 


2.1.2 Star formation and initial mass function 


The models differ in their treatment of star formation in galaxy 
discs (quiescent star formation). In the BAU model, the galactic star 
formation rate is calculated as the ratio of the total mass of cold gas 
to the characteristic star formation time scale, SFR = M g / T*, with 

'•■'■"GooEfO"'' 111 


where Vo is the circular velocity at the half-mass radius of the disc, 
T*o = 8Gyr and a* = —3. This prescription reproduces the ob¬ 
served gas fraction-luminosity relation at z = 0. 

In the FON model, the star formation time scale uses the pa- 
rameterisation 


raise ( Vo 3“* 

£* \200kms _l y 


( 2 ) 


where T^isc = cq/Vq is the disc dynamical time-scale and the free 
parameters e* = 0.0029 and a* = —1.5 were chosen mainly to re¬ 
produce the galaxy luminosity function at z = 0 ( {Bower et al.|2006) . 

In the LAG model, the star formation rate per unit area is as¬ 
sumed to be proportional to the surface density of molecular hy¬ 
drogen in the galactic disc. The amount of H 2 is obtained from the 
empirical relation of [Blitz & Rosolow sky 1 l |2006) which relates the 
fraction of molecular hydrogen to the pressure in the mid-plane of 
the galactic disc. 

The star formation prescription used in LAG has no free pa¬ 
rameters in the sense that they are calibrated against observations 
of star formation rates and gas surface densities in nearby galax¬ 
ies. The parameters are in effect measured and are only allowed to 
span a narrow range of values consistent with observational errors, 
narrower than the ranges within which other model parameters are 
typically allowed to vary. 

The stellar initial mass function (IMF) is assumed to be uni¬ 
versal in LAG and FON, of the form proposed by|Kennicutt|(|1983). 


2.1.3 Starbursts 

In all models, galaxy mergers can trigger bursts of star forma¬ 
tion, which transfer gas (and possibly stars) from the disc to the 
spheroidal component, where this gas participates in star forma¬ 
tion. The merger events are instantaneous in the models and lead 
to either an elliptical galaxy (for a merger of galaxies of similar 
masses) or to a corresponding increase in the galaxy mass if the 
merging satellite galaxy has a mass much smaller than the central 
galaxy. 

In LAG and FON, material is transferred to the spheroidal 
component, triggering a starburst also in the case when the galac¬ 
tic disc becomes dynamically unstable according to the bar stabil¬ 
ity criterion of |Efstathiou, Lake & Negroponte| {T982| l. The BAU 
model does not allow for any disc instabilities. 


2.1.4 Stellar and AGN feedback 

The models adopt similar parameterisations of the outflow due to 
supernovae and stellar winds, with a mass loss rate from the disc 
that depends on the star formation rate (a measure of the energy 
input into the ISM) and some power of the circular velocity of the 
disc (a measure of the depth of the potential well). The mass ex¬ 
pelled is reincorporated into the disc at a later time. In the BAU 
model, this happens once the dark matter halo has doubled its mass, 
whereas in LAG and in FON, this occurs after a time t = fdyn/Orh* 
where t& yn = ^Virial Avirial is the dynamical time of the halo and 

Or h = 1.26. 

On the other hand, the models differ significantly in their as¬ 
sumptions about the processes that regulate the bright end of the 
galactic luminosity function. In BAU, galactic super-winds are as¬ 
sumed to eject material from haloes that have low circular veloc¬ 
ity (with no subsequent reincorporation), at a rate proportional to 
the star formation rate. This reduces the baryon fraction in more 
massive haloes, which themselves are not directly subject to su¬ 
perwinds, thereby reducing the rate at which gas cools in these 
haloes. In LAG and FON, the growth of supermassive black holes 
in the centre of the galaxies is modelled and associated with an 
AGN feedback contribution. The cooling flow is interrupted once 
the black hole’s Eddington luminosity exceeds a multiple of the 
cooling luminosity. 


2.1.5 Ram-pressure stripping 

When a galaxy becomes a satellite galaxy (after the merger of dark 
matter haloes), it looses part or all of its hot gas halo to the central 
system by ram-pressure stripping. 

In BAU and LAG, this is modelled by transferring all of the 
hot gas halo of the satellite to the central galaxy, i.e., one assumes 
that the hot gas is completely and instantaneously stripped. In the 
FON, ram-pressure stripping is modelled more accurately, without 
assuming the instantaneous stripping of the gas. Instead, the pre¬ 
scription obtained from hydrodynamical simulations of |McCarthy| 
|et al.| ( [2008) is followed. 
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Table 1. Summary of the quantities used and notation. 


Type of parameter 

Notation Meaning 

Defined in or 
the value adopted 

Galactic properties obtained 

M g 

Cold gas mass of a galaxy 



from the semi-analytic models 

M* 

Stellar mass of a galaxy 



of galaxy formation 

**50.out 

Half-mass radius of the galactic disc 




To 

Circular velocity of the galactic disc at the half-mass radius 




SFR 

Star formation rate in the disc 



Quantities estimated in this paper 

h 

Average scale height of the galactic disc 

Eq.j 

0) 


P 

Average gas density in the galactic disc 

Eq. 

Z) 


ro 

Corrected half-mass disc radius 

Eq. 

a 


a 

Angular velocity of the disc 

Eq. 



s 

Maximum rotational shear 

Eq. 

m 


Vad 

Local outflow speed 

Eq. 

19) 

Adopted parameters 

lo 

Characteristic length scale of the turbulence 

0.1 kpc 


vo 

Root-mean-square gas velocity dispersion in the disc 

lOkms -1 


a 

Number of contributions to the interstellar pressure 

4 (Eq.|A6| 


Rk 

Ratio of turbulent diffusivities of the mean helicity and large-scale magnetic field 

0.3 (Eq^Tsj 

Computed quantities 

R„ 

Outflow magnetic Reynolds number 

Eq. 

27) 


D 

Dynamo number 

Eq. 

25)) 


D c 

Critical dynamo number 

Eq. 

29) 


B 

Steady-state large-scale magnetic fields strength 

Eq. 

30) 


b 

Steady-state random magnetic field strength 

Eq. {24) 



Figure 1. The late-type galaxy stellar mass function of the SAMs used. 
Different colours correspond to different SAMs as labelled. 


2.2 Late-type galaxies and their properties 

Magnetic fields ordered at the galactic scale only occur in late-type 
galaxies (since they have significant rotation), so we focus on these 
galaxies in what follows. We classify a galaxy as late-type if its 
bulge stellar mass accounts for less than half the total galaxy mass. 

Details of the galaxy formation models affect the distribution 
of galaxy stellar masses. In Fig. [T] this is shown for the explored 
models. 


2.2.1 Interstellar turbulence 


Interstellar gas is vigorously driven into a turbulent state by the 
injection of kinetic energy, mostly from supernova explosions (a 
general consequence of star formation), but also from stellar winds 
and other, less important sources {Mac Low & Klessen[|2004| 
|Elmegreen & Scalo||2004| |Scalo & Elmegreen 2004) . The root- 


mean-square (rms) turbulent speed is close to the sound speed in 
the warm gas (T ~ 10 4 — 10 5 K) at vo — 10kms~* and exhibits 
weak variations within and between spiral galaxies (Tamburro e t al.| 
|2009[ ). The turbulent scale Iq ~ 0.1 kpc is controlled by the size of 
a supernova remnant when its expansion velocity decreases to the 
sound speed in the ambient warm gas: the expanding supernova 
shell drives flows in its wider environment only from this stage on¬ 
wards. 


The small variability in the turbulent speed may seem surpris¬ 
ing as one might expect a dependence of vo on the star formation 
rate. Indeed, detections of such a dependence through HI observa¬ 
tions have been reported ( |Dib et al.|2006| ). However, both the in¬ 
terpretation of such observations and the connection between star 
formation rate and turbulence are far from straightforward. Apart 
from driving turbulent motions, supernova activity is responsible 
for the hot phase of the interstellar gas, and the filling factor of the 
hot phase increases with SFR ( |de Avillez & Breitschwerdt|2004| l. 
The energy released by supernovae is distributed among several 
such channels (including radiative cooling, acceleration of cosmic 
rays, etc.), and it is not obvious that an enhanced star formation rate 
would necessarily lead to a larger turbulent velocity. The fact that 
the turbulence observed in the warm interstellar gas is transonic is 
likely not to be a coincidence, but rather a result of non-linear feed¬ 
back between star formation and turbulence. If the turbulent driving 
becomes stronger because of an increase in the SFR, the turbulence 
initially becomes supersonic, leading to a rapid dissipation of ki¬ 
netic energy into heat in the shocks, so that the turbulent kinetic 
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energy is reduced, being converted into thermal energy of the warm 
gas (the speed of sound in the hot gas is of order 100kms~ 1 , so the 
warm gas absorbs most of the extra kinetic energy). As a result, the 
gas temperature increases until the turbulence becomes transonic 
and the system reaches a new (quasi-)steady state. However, the 
gas cools radiatively, on a cooling time scale (of a few 10 3 yr for 
temperatures near 10 4 K, assuming solar metallicity) that is shorter 
than the typical time scale of variations in the star formation rate 
(of order 10 6 yr). Thus, the warm gas can adjust itself to the varying 
star formation rate remaining at a nearly constant equilibrium tem¬ 
perature of T = 10 4 — 10 5 K which is known to be rather insensitive 
to the heating rate because of the efficient cooling. As a result, the 
extra energy supplied by the supernovae is more plausibly lost as 
radiation and, more importantly here, increases the fractional vol¬ 
ume of the hot gas. This enhances the associated outflow from the 
gas disc but does not increase the turbulent velocity. In other words, 
the turbulent velocity in the warm gas is controlled by its sound 
speed, i.e., by the balance of heating and cooling which, ultimately, 
is determined by its chemical composition (the metal abundance). 

Outflows (fountains or winds) driven by supernova activity en¬ 
train significant amounts of warm and cold gas (e.g., |Gent et al.| 
|2013[|McCourt et al.|2015) , so it is not surprising that observations 
at a relatively low resolution (1 kpc or larger) can pick up a broad 
range of velocities dominated by outflows rather than the turbulent 
motions in the warm gas. 

To conclude, we keep the turbulent velocity and scale fixed 
and independent of the star formation rate at vo = lOkms -1 and 
/o = 0.1 kpc. However, the outflow speed is a sensitive function of 
the star formation rate in our model (Section|2.2.5|l. 


2.2.2 Derived galactic quantities 

The calculation of the magnetic field in a galaxy that is described in 
Section[3]requires knowledge of the average gas density, p, and the 
average disc height, h. In this section we derive simple expressions 
for these two quantities as functions of galaxy properties output by 
SAMs (see Table [T] for a summary of the notation and quantities 
involved). 

Assuming that the surface density of the disc of the galaxy is 
well described by an exponential profile, the height of a gaseous 
disc in hydrostatic equilibrium is given by (see Appendix [A] for a 
derivation) 


h{r) 


a wl 


3 G(M* + M gas ) 


e. r / r » 


(3) 


where vo is the turbulent velocity of the disc gas, r s is the radial 
scale length of the disc and a is the number of pressure contribu¬ 
tions to the support of the gaseous disc. We use a = 4, assuming 
equal contributions from thermal, turbulent, cosmic-ray and mag¬ 
netic pressures. 

Equation 0 leads to a divergent galaxy volume unless it is 
truncated. We choose the half-mass radius as the galactic boundary, 
r o = i %r s with t, ~ 1.678, so that the volume of the galaxy follows 
as 


V = 2n [ h[r)r dr = 2nXrkh{Q ), 

Jo 


(4) 


where 


X = 


JO 


xs^ x dx 


1.645. 


h = 


The average disc scale height can be defined as 

V 


av o r o 


2 Ttfij 


3 G(M* + M gas 




which yields 


h ~287 pc 


ro 


v 0 


(lOkms^ 1 

l 


10 kpc 

a / M*+M gas 

X 4 V lO n M 0 
The average gas density is then obtained as 


P = 


Mi 


gas 

2V 


GM gds (M gas T A/* 


3^ 2 


CCVnr. 


,2 A 


o'o 


4kX 


(5) 

(6) 


(7) 


or 


p ~1.06 x 10 23 gem 3 



f A/gas \ / A/gas + Af* \ / Vo 

v 1 O 1 O M 0 J V 1 O U M 0 J V 10 kms-> 


-2 


(8) 


The factor 1/2 in equation 0 allows for the fact that ro in 
equation 0 is the half-mass radius, thus enclosing half of the total 
mass. 


2.2.3 Disc size 


Reproducing the observed sizes of galactic discs is a long stand¬ 
ing problem in semi-analytic models of galaxy formation (e.g. 
|Gonzalez et al.||2009} . Because of the strong dependence of the 
derived quantities on the disc size, we have renormalised the disc 
half-mass radii so that the medians of the predicted half-mass radii 
match the fit obtained by [Dutton et al.| ( [20 1 1) > for the observed re¬ 
lation between the half-mass radius and the stellar mass of SDSS 
galaxies at z = 0, 


GO.ernp ( 2 — 0) — r/) 



M*V 

M D J 


(Cl-Cl)/C3 


(9) 


where 


(ci,C 2 ,log 10 Af D ,log I0 ?- o ,C 3 ) = (0.2,0.46,10.39,0.75,1.95). 


Thus, the corrected half mass radius of a disc galaxy in a par¬ 
ticular mass bin is related to the original disc size in the output of 
the SAM, /"5o,out through 


r o(z) 


GO.(ml (2) 


GO.cmp 0) 

G0,median (o 0) 


( 10 ) 


This procedure corrects the median values of the disc radii in 
each mass bin, preserving the dispersion of disc sizes predicted by 
the models and the redshift evolution of the mass-size relation. In 
Fig. 0 we show the original and corrected size-mass relations for 
each SAM. 


2.2.4 Rotation and velocity shear 

An estimate of the disc angular velocity, £2, is necessary in the mod¬ 
elling of galactic magnetic fields. This is computed at the half-mass 
radius, ro, using 

Vo 

a = Q.o = —, ( 11 ) 

ro 

where Vq is the disc’s circular velocity at r = tq. 
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FON 



LAG 



Figure 2. Disc size-stellar mass relation for the three SAMs used here. The red curves and corresponding shaded areas show to medians and the range 
between the 15th and 85th percentiles for the raw model output. The blue curves and corresponding shaded areas show the same data after the rescaling to 
make the medians match the fit to observational data given by equation 0- Each panel shows a different model as labelled. 


The next quantity to compute is the shear, 5 = rdSl/dr. We 
employ the maximum value of S in our calculations. It can be 
shown that the maximum shear in a purely exponential disc is re¬ 
lated to the angular velocity through 

5 sa —0.76f2o • (12) 


that we have assumed that the OB associations can be treated as 
identical (i.e., they share a common equivalent mechanical lumi¬ 
nosity and the SN rate within them is approximately the same). 

The interesting quantity, however, is the average of v s h over 
the whole galaxy, v s h, taking into account the rate of occurrence of 
OB associations. This can be estimated as follows 


2.2.5 Galactic outflows 


The mean-field dynamo, responsible for the generation of a large- 
scale galactic magnetic field, relies on the removal of small-scale 
fields (and their magnetic helicity) away from the dynamo region 
in order to saturate in a steady state where the field has a strength 
comparable to that observed, of several microgauss (Brandenburg 
|& Subramanian|20 05 ). As discussed by Shukur ov et al.H2006) , in 
galaxies this can be achieved by the advection of magnetic fields 
by the outflow of the hot gas (and its entrained colder components) 
from the disc, as well as by turbulent diffusion i jKleeorin et al.|2000] 
and references therein). 

While the SAMs have internal prescriptions to compute the 
amount of gas which outflows from the disc due to SNe powered 
galactic winds, the gas carried by these winds is thought to be re¬ 
moved completely from the galaxy, and becomes available for rein¬ 
corporation into the hot halo only after time scales larger than dy¬ 
namical time of the dark matter halo. These strong winds, however, 
correspond only to a fraction of the total outflow. 

The outflow relevant for the magnetic field evolution is the 
one associated with galactic fountains, which leads to the removal 
of hot magnetised gas from the mid-plane of the disc but with a 
possibly very short re-accretion time-scale. This is not considered 
explicitly in the SAMs. We model this process as follows (see also 
|Lagos et al. 12013) . 

Supernovae tend to cluster in OB associations, regions of ac¬ 
tive star formation. The large, energetic bubbles of hot gas (super¬ 
bubbles) produced by dozens of supemovae and stellar winds in 
an OB association can more readily break out from the galactic 
disc. In the framework of the superbubble model of |Mac Low &| 
|McCray| ( j 1988| ), the speed of the shock front at the top of an ex¬ 
panding superbubble at it breaks-out is given by (see Appendix [b] 
for derivation) 

- 2/3 


v s h = 4kms 


”0 


- 1/3 


h 


(13) 


. 1 cm 3 J 1 kpc 

where iiq is the mean number density of the gas in the disc. Note 


l ; sh 


= v s h 


^OB 

Agal 


(14) 


where A ga ] = 71 ? g is the surface area of the galaxy. The total area 
filled with OB associations is 


Aqb = NoBtt(2h) 2 , 


(15) 


since the radius of an OB association at the break out is about 2h 
( |Mac Low & McCray|1988) . The number of OB associations can 
be obtained from the frequency of SNe occurrence in OB associa¬ 
tions, Vsn,ob. and the rate of supernovae occurrence in a single OB 

association, Vsniqb i 


Nob 


Vsn,ob 

vsn.iob 


(16) 


The frequency of SN occurrence in OB associations is ap¬ 
proximately 70 per cent of the total SN frequency |Tenorio-Tagle &| 
|Bodenheimer| j 1 988}, i.e. Vsn.ob = /ob v sn with /qb ~ 0.7. The 
overall SNe frequency relates to the star formation rate through 
Vsn = i?snxSFR. where T]sn = 9.4 x 10~ 3 Mq for Kennicutt’s 
IMF l |Kennicutt[ 1983) . 

The frequency of supernovae within a single OB associa¬ 
tion can be found from the number of SNe that typically oc¬ 
cur in one OB association, Asn,iob ~ 40 (Heiles|1987|, and the 
typical lifetime of an OB association is tps = 3 x 10 6 yr, i.e., 
v SN.lOB— toB^SN.lOB- Thus, we find that the number of OB asso¬ 
ciations in a galaxy at a given epoch is related to its star formation 
rate through 


Nob 


490 


SFR \ 
M 0 yr-V 


x 



tQB 

10- 3 Gyr 


Nsn.IOB 

40 



-l 


(17) 
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Using equations (| 1 3|) (| 17[> one obtains 



This is the gas speed in (and immediately behind) the shock front 
of a supperbubble. At the time scales involved, the magnetic field 
can be assumed to be frozen into the gas and thus is lost from the 
disc together with the gas. Therefore, the quantity of interest is the 
mass-averaged speed, v a( j, such that the surface density of mass loss 
from the disc is pv a( j: 

v ad — Uh — i (19) 

P 

where p/, ta 1.7 x 10~ 27 gcm~ 3 is the density of the hot gas and p 
is the average interstellar gas density. 


3 MODELS OF GALACTIC MAGNETIC FIELDS 


A partially ordered magnetic field B in the turbulent interstellar gas 
can be conveniently represented as the sum of a large-scale, B, and 
a fluctuating, b, components, 

B = B + b, b = 0, 


where a bar denotes the ensemble or any other suitable average (see 
|Gent et al.|2013[ for a discussion of averaging procedures). 

The growth time of the large-scale magnetic field due to the 
mean-field turbulent dynamo jMoffatt|1978||Ruzmaikin et al.| 1988} 
can be estimated as jji et al.|2014fr 

IMF — 2 —(D c —Df 1 / 2 , \D\ > |Z> C |, (20) 

where the magnetic diffusivity p, dominated by the turbulent con¬ 
tribution, is estimated from mixing-length theory as 

P — jhvo, 


with Zq and vq the turbulent scale and speed, respectively, the dy¬ 
namo number D (that quantifies the strength of the dynamo action) 
can be written using equations © and \26\ as 



with 5 being the rotational shear rate, and D c is a critical value 
of the dynamo number such that magnetic field can be maintained 
only if D < D c (D < 0 so long as £2 decreases with r). We note that 
h 1 /r\ is the turbulent diffusion time across a layer of a scale height 
h. Equation ( |20| l is quite accurate for the range of | D of interest in 
applications to spiral galaxies (see the discussion in |Ji et al.|2014) . 
For a flat rotation curve, S = —£2, and \D\ |D C |, which is true in 

the inner parts of most galaxies, we obtain D ~ — 10(Z?£2/vo) 2 and 


2 h 

™ f ^Zo£2 

" 4xl ° 8yr ((T5l^)( 25k^s/kpc ) 


to 

0.1 kpc 


( 21 ) 


( 22 ) 


normalised to the parameter values corresponding to the Solar 
vicinity of the Milky Way. In the inner parts of slowly rotating 
galaxies and in outer regions of normal disc galaxies, Tmf is shorter 
than the galactic evolution time scale T e ~ 10 9 yr. 


Random magnetic fields grow even faster due to fluctuation 
dynamo action (note that this dynamo action is distinct from the 
more commonly discussed mean-field dynamo), on a time scale 
shorter than the turnover time of the largest turbulent eddy, To ~ 
Zo/vo ~ 10 7 yr. 

The exponential growth on the time scales specified above 
halts (the dynamo action is said to saturate) as soon as the Lorentz 
force becomes comparable to other forces on the relevant length 
scale, of order 1 kpc, for the large-scale magnetic field and a frac¬ 
tion of the turbulent scale for the random field ( |Bhat & Subrama-] 
nian| |2013| l. After that, magnetic fields remains in a statistically 
steady state. The magnitudes of B and b in this state are discussed 
further in this section. A recent review and a suite of formulae de¬ 
scribing nonlinear mean-field dynamos can be found injChamandy] 
|etak] (2014). 

Thus, as a first approximation, we can assume that: (i) galac¬ 
tic magnetic fields adjust themselves instantaneously to the evolv¬ 
ing galactic environment (i.e, Tmf,To <C T e ); and (ii) they always 
remain in a statistically steady state of a saturated dynamo. These 
assumptions are obviously crude. However, they are sufficient for 
a first exploration of the effects of magnetic fields on galaxy evolu¬ 
tion if the exploration then continues to include the effects of finite 
dynamo time scales. The latter will be our subject elsewhere. 


3.1 Random magnetic fields in spiral galaxies 

Any random flow of electrically conducting fluid is a dynamo (i.e., 
it amplifies a seed magnetic field exponentially fast) provided the 
magnetic Reynolds number due to Ohmic magnetic diffusivity, 
Rm = Zovo/p, exceeds a certain critical value, R m , C r of order 100 in 
an incompressible flow. This type of a dynamo is known as the fluc¬ 
tuation dynamo ( jZeldovich et al.| l990j [Schekochihin et al.j2002[ 
|Brandenburg & Subramanian[2005[l. 

The growth time of the random magnetic field in a vortical 
random velocity field of a scale Z is as short as the kinematic time 
at that scale, T ~ Z/v(Z) with v(Z) the rms random speed at the scale 
Z. In a turbulent flow with a sufficiently shallow power spectrum, T 
is shorter on smaller scales. For example, T « Z 2 / 3 in a flow with 
the Kolmogorov spectrum v(Z) « Z 1 / 3 . Therefore, the magnetic en¬ 
ergy spectrum peaks at small scales during the exponential growth 
(kinematic) stage, and then plausibly settles to a form similar to 
that of the kinetic energy spectrum in the statistically steady state 
(a saturated, nonlinear dynamo). The eddy turnover time in spi¬ 
ral galaxies, where Iq ~ O.lkpc and i’o — lOkms -1 , is as short as 
To — Iq/vq ~ 10 7 yr even at the energy scale of the interstellar tur¬ 
bulence. Thus, the fluctuation dynamo can rapidly produce random 
magnetic fields in the ISM i jShukurov|2007} . 

Compressibility hinders the fluctuation dynamo, reducing the 
growth rate of the random magnetic field. In the extreme case of 
sound-wave turbulence, the longitudinal nature of the fluid motions 
reduces the probability of three-dimensional twisting and folding 
of magnetic lines, an essential element of the amplification mech¬ 
anism (as in Zeldovich’s stretch-twist-fold dynamos -|Zeldovich,| 
|Ruzmaikin & Sokoloff|1983||1990| >. As a result, the growth time of 
the rms magnetic field is of order T ~ _-#~ 4 Zo/vo for 1, where 

is the Mach number (Kazantsev et al. 1985). A generic com¬ 
pressible turbulence inherits this feature, producing a higher thresh¬ 
old value of R m and slower growth of magnetic energy (Haugen 
|et al,|2004l|Federrath et al.|2011||Gent et al,|2013) . We note, how¬ 
ever, that the velocity field in such flows is a mixture of solenoidal 
and compressible parts ( |Moss & Shukurov) 1 996) , with their rela¬ 
tive contributions depending on the Reynolds number, Mach num- 
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ber, numerical resolution, etc. Haugen et al. 2004, Mee & Bran¬ 
denburg 20061. Haug en~et al.| ' 2004| > find, in their simulations of 
dynamo action in a compressible random flow, that R mcr increases 
from about 35 at = 0 to 60 at = 1 and 80 at = 2.5 (for the 
magnetic Prandtl number P m = R m /Re = 1 with Re the Reynolds 
number) but suggest that R mjC r may vary little as the Mach number 
increases further, especially for a_large magnetic Prandtl number 
(as in the interstellar gas). 


Federrath et al. 


(20111 find roc „#-V 3 


for 10 < < 20 in their simulations of isothermal compressible 

random flows. 

However, slower growth and a higher dynamo threshold in a 
compressible random flow, as compared with incompressible flows, 
may not in fact represent a practical problem since the value of the 
magnetic Reynolds number is usually very high in astrophysical 
plasmas, and the turbulent kinematic time scale Iq/vq is very short 
in comparison with any global time scale in galaxies. What is more 
important is the steady-state magnitude of the magnetic energy den¬ 
sity relative to the kinetic energy density of turbulence. According 
to |Federrath et al.] < |201 1| ), isothermal compressible random flows 
produce lower magnetic energy density /8 n relative to 

the kinetic energy density E± = l,pv g, with E m /E^ ~ 10~ 2 -10~ 3 
at 2 < < 20 than a purely solenoidal velocity field. Thus, gas 

compressibility is detrimental to the fluctuation dynamo in both the 
kinematic and saturated regimes. 

The magnetic field produced by the fluctuation dynamo is in¬ 
termittent, and can be described as a statistical ensemble of mag¬ 
netic flux ropes and sheets whose coherence size (length or ra¬ 
dius of curvature) is of the order of the flow correlation length, 
Iq, whereas the rope thickness is, in t he kinematic dynamo, of the 


order of the resistive scale, IqR, 


- 1/2 


i Zeldovich et al. 


19901. 


Wilkin 


|et al.[ < [2007) l show that flux ropes become progressively dominant 
as R m increases in the kinematic regime. 

Kinematic fluctuation dynamos are well understood, but the 
nonlinear, statistically-steady state remains, to some extent, contro¬ 
versial. Simulations of fluctuation dynamos in driven random flows 
suggest that magnetic energy density within the ropes is close to 
equipartition with the kinetic energy density ( jBrandenburg & Sub-| 
|ramanian|2005}, 


Seq = V Q {ATip) 1 / 2 . 


(23) 


A widely accepted model for the saturation of the fluctuation dy¬ 
namo by Subramanian (1999) (see also Subramanian et al.(2006), 
suggests that, in the steady state, the magnetic Reynolds number 
is renormalised to its critical value, so that the thickness of mag- 
netic ropes is estimated as d ~ h)R m ,a — 0.1 Zo for R m ,ci = 100 
independently of R m (see, however, |Schekochihin et al.|2002| ). For 
modest effective magnetic Reynolds numbers, R m ~ Rm, cr> mag¬ 
netic sheets may be predominant. Then the volume filling factor f B 
of the magnetic structures can be estimated assuming that there is 
one magnetic sheet per correlation cell of the turbulent flow: 

f B = j- — Rm,a — 0 . 1 , 

so that the root-mean-square magnetic field follows as 
b~f B B eq =f B v 0 (47zp) l F. (24) 


However, magnetic field outside the dominant magnetic structures 
can still contribute significantly to the observables, such as the ran¬ 
dom Faraday rotation measure i jBhat & Subramanian|2013| , so that 
this estimate should be applied judiciously. 


3.2 Regular magnetic fields in spiral galaxies 


The generation of large-scale magnetic fields in galaxies is de¬ 
scribed by galactic dynamo theory using the concept of a mean- 
field dynamo (see Be ck et al.|199 6, Shukurov 2007] for a review). 
Dynamo action in a rotating, stratified galactic gas layer is pro¬ 
duced by the joint action of the helical turbulent motion (via the 
so-called a-effect) and the differential rotation of the galactic disc. 
These two effects are quantified, respectively, by two dimension¬ 
less parameters, 


Ra — 


0()/j 

V 


and 


R<o — 


Sh 2 

1) 


The widely used aco-dynamo approximations applies where 
|Rid| Ra- In this approximation, dynamo action is controlled 
solely by their product, known as the dynamo number, 


D — RaRco ■ 


(25) 


The magnitude of the cc-effect can be obtained from the order-of- 
magnitude estimate 



(26) 


It is not quite clear how the mean-field dynamo enters a non¬ 
linear, steady state (for a review, see|Brandenburg & Subramanian| 
2005). In the most detailed and physically motivated theory avail¬ 
able, the growth of the large-scale magnetic field is limited by the 
conservation of magnetic helicity, so its steady-state strength de¬ 
pends on the rate at which the mean magnetic helicity of the random 
magnetic field is removed from the localisation region of the large- 
scale magnetic field (Del Sordo et al. 20131. In galaxies, this can be 
accomplished by galactic winds and fountain flows ( Shukurov et al. 
|2006|[Sur et al.|20 07 1. In this case an additional dimensionless pa¬ 
rameter enters the picture, quantifying the advection of magnetic 
helicity out of the galactic disc 



(27) 


Another mechanism that contributes to the partial removal of 
the small-scale magnetic helicity is its turbulent diffusion l jKleeorin| 
|et al.|2000||2002||2003| l. Allowing for this effect introduces another 
dimensionless number, the diffusivity ratio of the mean current he¬ 
licity, fc, and the mean magnetic field, rj: 

Rk = - ■ (28) 

V 

|Mitra et akj ( |2010| obtained R K ~ 0.3 independent of the Reynolds 
number. 


To summarise this model (see |Chamandy et al.|2014| for de¬ 
tails), a disc galaxy has an active mean-field dynamo (i.e., it is ca¬ 
pable of amplifying and sustaining a large-scale magnetic field) if 
the the magnitude of the dynamo number, equation ( |25| >, exceeds 
that of its critical value, 



(29) 


For \D\ < |D C |, the large-scale magnetic field decays on a time scale 
of order h 2 /r\. Otherwise, B is exponentially amplified up to the 
steady-state strength 

(£-i)(^ + ^), (30) 


where the magnetic pitch angle, defined as p = arctan [BjB^ 
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Figure 3. The predicted magnetic strength function (MSF), equation of the large-scale magnetic fields of spiral galaxies at z = 0 in different SAMs, 
as indicated by the key. Each panel corresponds to a different galactic stellar mass interval, indicated in the top left of each frame. The points correspond to 
binned model results. 


in terms of the cylindrical magnetic field components, can be ex¬ 
pressed in term of the previous quantities as 


1 

p = arctan ( — 
, Bm 


2D C 


7T 1 / 2 


with B eq defined in equation |23j, and 
C(p) = ( i-^cos 2 p 


(31) 


(32) 


4 RESULTS AND DISCUSSION 


In this section we identify and interpret robust features of galactic 
magnetic fields obtained from the galaxy formation models. Our re¬ 
sults refer to statistical properties of large-scale and random mag¬ 
netic fields in large samples of galaxies. We found that magnetic 
properties are different in galaxies of different masses and our re¬ 
sults are presented for characteristic galactic mass ranges, where 
we have calculated the mass using only the stellar content. 

The statistical properties of galactic magnetic fields are con¬ 
venient to describe in terms of the number density of galaxies with 
a particular magnetic field strength per logarithmic interval of the 
filed strength; in the case of the large-scale magnetic field, this vari¬ 
able, referred to as the magnetic strength function (MSF), is defined 
as 


W (B) 


dn 

d logS 


(33) 


4.1 The large-scale magnetic field 

Figure[3]shows the distribution of magnetic field strength, the MSF, 
for the large-scale field B in the local Universe (z = 0), each panel 
displaying a different galactic mass interval. The form of the MSF 
is sensitive to the galaxy formation model and, in each model, to 
the mass interval. 

The typical magnetic field strengths are similar in the FON 
and BAU models, 0 < log 10 (B/l/iG) < 1. Except for the highest 
mass interval, the MSF of the LAG model is clearly bimodal, with 
the first peak in the interval 0.2 — 0.8 pG and the second close to 
the what is predicted by the other models. 

Lower mass galaxies can have the strongest large-scale mag¬ 
netic fields. Strong large-scale magnetic fields are increasingly 
suppressed with increasing mass and for galaxy masses M* > 
IqIO.25 Mq there is a negligible number density of galaxies with 
large-scale magnetic fields stronger than ~ 10/lG. The reason for 
the suppression, discussed in more detail in section |4~4| is the in¬ 
crease in the outflow speed in these galaxies, associated with their 
higher star formation rates. 

Another factor that affects Fig.[3]there is the overall decrease 
in the number density of galaxies with the galactic mass. This is 
expected from the shape of the galaxy stellar mass function of late- 
type galaxies, shown in Fig. |T| However, as it will be discussed 
later, the decrease is also intensified by the relatively small fraction 
of massive galaxies with active dynamos. 

In Fig. [4] we show the fraction of galaxies with active large- 
scale dynamos (i.e., \D\ > |Z) C |) as a function of galaxy mass at z = 
0. The fraction of active dynamos in galaxies decreases with mass 
and all models predict that fewer than 40 per cent of galaxies have 
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Figure 4. The fraction of galaxies with active dynamos (i.e. \D\ > D c ) at 
Z = 0 as a function of galaxy mass in each model, as shown in the legend. 


active dynamos for M > 1O 1O M 0 . Fig. [^provides a more detailed 
picture where the fraction of galaxies with a magnetic field strength 
exceeding 0.1,1 and 10 /r G is shown as a function of mass for each 
galaxy formation model separately. 

The bimodal nature of the MSF of the LAG model in Fig. [3] 
is clarified by Fig. [6] where the MSF is shown separately for the 
central galaxies and their satellites: both the central galaxies and the 
satellites have unimodal MSFs but with maxima at different values 
of B in the LAG model. This is a consequence of the differences in 
the gas content of the satellite and central galaxy populations (see 
section|4~4|for details). 


4.2 The total magnetic field 

Fig. [7] shows the MSF of the total magnetic field comprising both 
the large-scale and random parts, B = ( B 2 + £> 2 )'/ 2 ; the MSF of 
the central galaxies is also shown with a dotted curve. The overall 
features of the large-scale magnetic field distributions of Fig.[3]can 
be seen in the total field as well. The most important differences 
occur at the two highest mass intervals where the MSF of the total 
magnetic field is significantly broader than the MSF of B. At these 
masses, the small-scale magnetic fields dominate. This is consis¬ 
tent with the small fractions of active large-scale dynamos in these 
galaxies, as shown in Fig. [5] 

In the LAG model, the tail of large-scale magnetic fields 
B < 0.1 flG is concealed by stronger random magnetic fields, pro¬ 
ducing a cut-off at a mass-dependent minimum field strength. On 
the other hand, in the FON model, the opposite happens: there is 
a population galaxies (both satellites and central) with negligible 
large-scale magnetic field which contains random magnetic fields 
b < 0.2 flG. Such a change does not occur in the BAU model where 
the total magnetic field has a distribution not dissimilar to that of 
the large-scale field alone. 


4.3 Evolution of galactic magnetic fields over cosmological 
time scales 

Fig. m shows the redshift evolution of the median value of the 
steady-state large-scale magnetic field in each galactic mass inter¬ 
val, as well as the 15th and 85th percentiles. Thus, at a fixed mass 
interval, the median of the large-scale magnetic field increases 
quasi-exponentially with redshift. The main reason for this increase 




log(MJM s ) 

Figure 5. The fraction of galaxies with B > 0.1, 1 and 1 0/iG at z — 0 as a 
function of the galaxy mass in each model, as specified in the legend. 


are the larger interstellar gas densities that occur at higher redshifts. 
One must, however, bear in mind that Fig.[8]reflects changes in both 
magnetic field strength and galaxy mass with redshift. 

The large-scale mag netic fields are expected to evolve on a 
time scale of 10 8 -10 9 yr i Beck et al7|| 1994 Shukurov|2007| . On 
the one hand, the dynamo time scale is shorter than the galactic 
evolution time, which allows us to use the steady-state strength of 
magnetic field in our estimates presuming that the dynamo action 
generates the large-scale magnetic field in a young galaxy and then 
adjusts itself instantaneously to the evolving galactic environment. 
This approximation is evidently acceptable in the inner parts of a 
galaxy where the rotational velocity shear is high and hence the 
dynamo time scale is short. However, the effects of a finite dynamo 
time scale can be significant in the outer parts of galaxies. 

To illustrate these arguments, consider again the time scale 
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Figure 6. The distribution of large-scale magnetic field strengths at z — 0; each panel shows different galaxy stellar mass interval, as labelled. Dotted curves 
show the MSF associated with satellite galaxies and continuous lines show the MSF of central galaxies. 




B [/xG] B [/xG] 




B [fiG] B [/xG] 

/_ 2 \ 1/2 

Figure 7. The distribution of the total magnetic field strength, B = ( B +b 2 ) , in different galaxy stellar mass intervals at z = 0, as labelled in each panel. 

Dotted lines show the corresponding distributions for the central galaxies alone. Different colours are used to indicate different SAMs as specified in the 
legend. 
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Figure 8. Redshift evolution of the large-scale magnetic field. The solid lines show the median field values. Dashed and dotted lines correspond to the 15 and 
85 percentiles, respectively. Different columns correspond to different models (as indicated in the top row) and different rows show the results for different 
galaxy mass intervals (as indicated). 


of the mean-field dynamo given by equation (20). Assuming, for 
the sake of illustration and by analogy with the Milky Way, that 
h = 0.2kpc at r = lkpc and £1 = VoA with Vy = 200kms -1 , we 
obtain for the inner galaxy 

™^ 2 ° Myr (lltpc) ’ 

having adopted D c m — 8 and D ~ 30 D c at r ~ 1 kpc. Thus, 
the e-folding time of the large-scale magnetic field can be as short 
as 20Myr in the inner galaxy. With h = 0.5kpc at r = lOkpc, and 
adopting D c — D = 4 for illustration, we have Tmf — 0.7 Gyr in 
the outer galaxy. The effective large-scale seed magnetic field pro¬ 
duced by the fluctuation dynamo is estimated by |Beck et a 1. | ( | 1994) 
as B s = 10~ 3 jUG; then the time required to amplify it to a strength 
6=1 flG is about 7TMF- These estimates suggest that the inner 
few kiloparsecs of galactic discs would host /xG-strong large-scale 
magnetic fields very soon after their formation. A galaxy formed at 
a redshift z = 10 would have its central kiloparsec magnetized with 
a large-scale magnetic field by z ~ 8, whereas a microgauss-strong 
magnetic field can build up via the local dynamo action by z ~ 1 at 
r = lOkpc (see also |Beck et al.|1994) . 

The above estimate is rather conservative since, due to the 
stronger dynamo action in the inner galaxy, the outer radius of 


the magnetized region increases approximately linearly in time 
at a speed Vf ~ iMoss et al. 1998); this expansion can be 
described as the propagation of a magnetic front driven by the 
mean-field dynamo action. Taking h = 0.3 kpc as a representa¬ 
tive value for the whole disc and r) = 10 26 cm 2 s _1 , we obtain 
Vf ~ 6(r/lkpc) _1 / 2 kms _1 , and the front reaches r = lOkpc in 
t = /q dr/Vf ~ 3.5Gyr. Thus, a magnetic front propagating out¬ 
wards from the inner galaxy can produce a strong magnetic field in 
the outer galaxy by z — 1.5, sooner than suggested by the above 
estimate of the local growth time at r = lOkpc. Effectively, the 
propagating magnetic front provides a strong seed magnetic field 
for the local dynamo action thus leading to a faster build-up of the 
large-scale magnetic field. 

We stress that turbulent magnetic fields can be produced on 
much shorter time scales anywhere in the galaxy (Section[3.1|l. 


These estimates of the redshift at which a spiral galaxy can de¬ 
velop a /tG-strong large-scale magnetic field are somewhat higher 
than those of Arshakian et al. (2009, 20111 since these authors did 
not consider the fact that the dynamo time scale is shorter in the 
inner galaxy: their estimates apply to the local dynamo action at 
r ~ lOkpc. 

Since the dynamo model used here does not include the dy¬ 
namo time scales explicitly, assuming that they are shorter than the 
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galactic evolution times, we only extend our results to z = 2.5. At 
this redshift, the outer radius of the part of the galactic disc occu¬ 
pied by 

the large-scale magnetic field of a microgauss strength is about 
6 kpc, allowing for local dynamo action alone. Extension of the 
model to include evolutionary equations for magnetic field will be 
published elsewhere. 

Figures |9| 1 1 0| and 1 11 1 show the MSF at different redshifts for 
each galaxy formation model. Except for a shift towards higher val¬ 
ues of B, there is not much variation in the shape of the MSF with 
redshift for the less massive galaxies. 

In the same figures, dotted curves show the MSF of central 
galaxies alone, suggesting that there is no significant change in 
the relative contributions of the central and satellite galaxies to the 
MSFs. 


4.4 Magnetic fields and the gas content of a galaxy 

The form of the MSF is closely related to the gas content of galax¬ 
ies. The gas density enters the reference magnetic field strength B eq 
through equation (|23j. Moreover, the efficiency of magnetic helic- 
ity advection, quantified by R„ of equation (|27j, depends on the star 
formation rate, which is controlled by the amount of cold gas in the 
galactic disc available to be converted into stars. 

Equation ( |30| l shows that galaxies with stronger star forma¬ 
tion, and hence larger R u , can have stronger large-scale magnetic 
fields. However, larger values of R„ also make the critical dynamo 
number higher, as shown by equation \29\ . As a result, there exists 
a star formation rate optimal for the mean-field dynamo action. In 
galaxies with large R u the dynamo action is suppressed. Because of 
their higher star formation rates, this affects massive galaxies es¬ 
pecially strongly, systematically reducing, with the galactic mass, 
both the strength of their magnetic fields (as can be seen in the 
plots of the MSF) and the fraction of galaxies with active dynamos 
as shown in Fig. [4] 


4.5 The importance of the ram-pressure stripping 

An important consequence of the non-linear relationship between 
the gas content, star formation and the magnetic field strength is 
the difference between the MSFs of central and satellite galaxies 
shown in Fig. [ 6 ] Because of the loss of hot gas by the satellites due 
to ram-pressure stripping, both the gas density and star formation 
rate in them decrease leading to weaker outflows, thereby lowering 
B. The effect is more evident for the FAG model because of the 
assumption of instantaneous ram-pressure stripping of the hot gas 
halo of the satellite (the so-called satellite starvation). For the FON 
model, the separation between the MSFs of satellite and central 
galaxies is smaller, reflecting the gradual ram-pressure stripping 
adopted in this model. 

In the BAU model, despite the same efficiency of the ram- 
pressure stripping as in the LAG model, the MSF has a maximum 
at nearly the same positions for both the satellite and central galaxy 
populations. The reason is the star formation law adopted in this 
model, which leads to a higher star formation rate in galaxies with 
smaller circular velocity, as it follows from equation[T] The increase 
in star formation rate in lower-mass galaxies leads to an increase in 
R„, compensating the effects of starvation. 


4.6 Magnetic field strength and star formation rate 


The dependence of magnetic field strength on the SFR is due to an 
explicit dependence of the large-scale field strength ( ]30) on the out¬ 
flow Reynolds number defined in equation \21\ , with v ac j given 
in equations ( | 1 9| > and < | 1 8| ). In addition, the SFR is involved in the 
mean-field dynamo action through the dependence of the critical 
dynamo number on R u in equation (29}. With the galactic outflow 
model used here (AppendixjB}, the outflow speed is proportional to 
the energy supply rate from supemovae in OB associations, 
and thus to SFR 1 / 3 . However, we then assume that all OB associ¬ 
ations have the same energy supply rate, so that the outflow speed 
is proportional to the SFR. As a result, equation (30} predicts that 
B SFR 1 / 3 in galaxies with a strong outflow (where R„ > n 2 R K ), B 
is independent of the SFR if the outflow is weaker, and the dynamo 
action is suppressed completely by very strong outflows, such that 

l D l < |A:|. 

The overall dependence of the total and large-scale magnetic 
fields strengths on the SFR obtained in our model is shown in the 
bottom and top rows of Fig. [12] respectively. The large-scale 
magnetic field strength is roughly independent of the SFR and 
the model predictions are compatible with observations of nearby 
spiral galaxies compiled by |Van Eck et al.| < [2015| ). The second 
row shows the variation of total magnetic field strength with the 
SFR. The model predictions are compatible with observations of 
both nearby spiral galaxies and Local Group dwarf irregular galax¬ 
ies ( jChyzy et al.||20lT} . For comparison, we also show (the 
dashed grey curve) the frequently adopted power law dependence 
B °c SFR"' with m ~ 1/3. A power law dependence is compatible 
with the predictions of the BAU and LAG over narrow SFR ranges 
and of the FON model over the whole SFR range. 

When interpreting the lower row of Fig. [72] one ought to bear 
in mind that, since the total magnetic field is dominated by the ran¬ 
dom component and the turbulent velocity is fixed in our models 
(see Section [3T} , the emerging relation between the field strength 
and the SFR is a consequence of the non-linear relation between 
gas density and SFR and not of the dependence of the turbulent 
velocity on the SFR. 

|Schleicher & Beck| ( |2013|) argue that inverse Compton scatter¬ 
ing of cosmic ray electrons can dominate their synchrotron losses 
from redshift z ~ 3 upwards, thus affecting the observability of 
galactic magnetic fields at high redshift. We note, however, that 
these authors use interstellar gas parameters obtained by |Thomp2| 
[son et aljj (2005} for starburst galaxies with the surface density 
of star formation E* — lO 3 M 0 kpc _ 2 yr _I , where vertical support 
against gravity is provided by radiation pressure. However, they 
apply this model to galaxies with surface density of SFR of order 
E* = 0.1 Mq kpc ~ 3 yr ~ 1 where the properties of the interstellar gas 
are rather different. Perhaps more importantly, Schleicher & Beck| 
12013 ' employ a simple relation B oc E.J / 3 for the total magnetic 
field strength (dominated by the random magnetic field). As dis¬ 
cussed above, any relation of B to E* is likely to be more compli¬ 
cated and can definitely be different in galaxies of different masses 
and different star formation rates. As discussed in Section |2.2.1| the 
main reason for this is that the turbulent speed in the interstellar gas 
is self-regulated to be comparable to the sound speed (if the energy 
supply is high enough), because supersonic turbulence quickly dis¬ 
sipates the excess kinetic energy into heat thus increasing the speed 
of sound. In galaxies with strong star formation, a progressively 
larger part of the energy supplied by supernovae and stellar winds 
is channeled into galactic outflows rather than turbulence. Since the 
hot gas leaves the disc at time scales shorter than even the fluctua- 
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Figure 9. The evolution of the MSF in the case of the FON model. The dotted curves correspond the MSF of central galaxies. Different colours indicate the 
predictions for different redshifts as labelled. 
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Figure 10. The evolution of the MSF for the BAU model. The dotted curves correspond the MSF of central galaxies. Different colours indicate the predictions 
for different redshifts as labelled. 
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Figure 11. The evolution of the MSF for the LAG model. The dotted curves correspond the MSF of central galaxies. Different colours indicate the predictions 
for different redshifts as labelled. 


tion dynamo time, both large-scale and random magnetic fields of 
galactic discs are generated in the warm gas where the speed of 
sound is of order lOkms -1 . This strongly nonlinear effect seems 
to preclude any universal dependence of the turbulent speed and 
magnetic field strength on the SFR applicable beyond relatively 
narrow ranges of galactic mass and SFR. This makes the conclu¬ 
sions of |Schleicher & Beck| ( |2013| l questionable. The dependence 
of galactic magnetic field strength on star formation rate is usu¬ 
ally discussed in connection with the radio-far-infrared correlation 
(e.g. Yun et al. |2001| and references therein), where this depen¬ 
dence is just one element of a complex physical system. It is not 
our intention to discuss or explain that correlation, even though our 
conclusions are relevant to the discussions of its nature. 

Since galaxy mergers are instantaneous in SAMs, we cannot 
include any effects of the mergers on galactic magnetic fields. How¬ 
ever, the starburst produced by a merger is a part of the SAMs; 
both large-scale and turbulent magnetic fields are enhanced dur¬ 
ing a starburst in our model. The effects of galaxy mergers on the 
mean-field dynamo can be nontrivial (Drzazga et al. 2011; Geng 
|et al.|2 012 Kotarba et al. 2011), but their inclusion would require a 
more detailed model of both galaxy formation and dynamo action 
that those used here. 

4.7 Advection and diffusion of magnetic helicity 

Table[2]shows that most galaxies with an active mean-field dynamo 
have R u < n. Thus, the strength of the large-scale magnetic held, 
given in equation (30) , is dominated by the diffusion of magnetic 
helicity for the value of R K used, R K = 0.3. This reduces signif¬ 
icantly the effect of galactic outflows on the mean-held dynamo 
action. 


Table 2. Fraction of galaxies with active dynamos and R„ <n in the galaxy 
formation models explored. 


Mass interval 
logM*/MQ 

LAG 

percentage 

BAU 

percentage 

FON 

percentage 

8.00-8.75 

85 

88 

93 

8.75-9.50 

89 

89 

95 

9.50-10.25 

92 

87 

92 

10.25-11.00 

92 

89 

90 


Equation (30) has been derived by Chamandy et al. (2014) 
with the diffusive flux of the magnetic helicity a m (proportional to 
the mean current helicity of the random magnetic held, b ■ V x b) 
obtained assuming that the scale height of a m is equal to that of the 
ionised gas. Almost nothing is known about the spatial distribution 
of the mean current helicity even from numerical simulations, not to 
mention observations. The numerical simulations of IGressel et al.l 
(2008) (see also |Gressel et al.|2013) suggest that the scale of the to¬ 
tal a-effect across the gas layer can be signihcantly larger (at about 
lkpc in the Solar vicinity of the Milky Way) than the gas scale 
height, perhaps due to the contribution of a m . If this is the case, the 
value of R K is strongly overestimated. Then the dynamo model un¬ 
derestimates the role of galactic outflows and, hence, star formation 
on the large-scale magnetic held. The relative roles of galactic out¬ 
flows and diffusion in the nonlinear mean-held dynamos requires 
further careful study under realistic conditions for spiral galaxies. 
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Figure 12. The dependence of the large-scale field strength, B (upper row), and total magnetic (lower row), B, on the star formation rate (SFR) predicted by 
each model. The solid lines show the median magnetic field values while the shaded areas show the range between the 15th and 85th percentiles. The square 
data points and triangles correspond, respectively, to the observed values and upper limits obtained by jChyzy et al.|j2011) for dwarf irregular galaxies in the 
Local Group. The diamond-shaped data points were calculated from the compilation made by jVan Eck et al.|j2015| their Table 2): when the surface density 
of star formation, £sfr was absent, we adopted the value corresponding to the innermost galactic radius; the total SFR was then obtained by integrating the 
surface density of SFR, Zsfr, and the overall magnetic field was taken as the area-weighted average of the reported values. The dashed grey curve in the 
bottom row shows the relation B SI R 1 ' with an arbitrary choice of the intercept (see text). 


5 SUMMARY AND CONCLUSIONS 

By coupling a nonlinear mean-field dynamo model with three well- 
established semi-analytic models of galaxy formation, we have de¬ 
veloped a framework to predict the strength of large-scale (global) 
and small-scale (turbulent) magnetic fields in evolving spiral galax¬ 
ies from the time when their bulge to disc mass ratio reduces be¬ 
low 1/2 and they are classified as late types. We present our results 
for various galactic stellar mass ranges selected according to the 
physical properties of the galaxies and we consider the evolution 
of galactic magnetic fields with redshift. Our main assumption is 
that the steady-state strength of the magnetic field is established 
instantaneously as the host galaxy evolves. This assumption is not 
questionable for the turbulent magnetic fields as the time scale of 
the fluctuation dynamo producing them can be as short as lOMyr. 
The time scale of the large-scale (mean-field) dynamo is longer at 
about a few Gyr in the main part of a spiral galaxy but can be two 
orders of magnitude shorter in the inner galaxy. This time scale is 
at least marginally shorter than the galactic evolution time, which 
justifies our assumption. 

We find that the choice of the galaxy formation model strongly 
affects the number density of galaxies that host a large-scale mag¬ 
netic field of a given strength (called the MSF in the text) in most 
mass ranges. In other words, the probability distribution of the 
strength of the large-scale magnetic field in a representative sample 
of galaxies is sensitive to the galaxy formation model. We discuss 


how the shape of the MSF is related to the physical processes af¬ 
fecting the dynamo action and how these depend on the galaxy for¬ 
mation model. In particular, the ram-pressure stripping of the gas in 
the haloes of satellite galaxies results in quite different typical field 
strengths in satellites and their central galaxies, which can lead to a 
pronounced bi-modality of the MSF. Our experience thus demon¬ 
strates the possibility of using observations of galactic magnetic 
fields to constrain galaxy formation models. 

Our results are also relevant to planning observational studies 
of galactic magnetism at high redshifts. Galaxies that host a large- 
scale magnetic field can produce a signal in polarised synchrotron 
emission, even when they are unresolved: the degree of polarisa¬ 
tion depends on the orientation of the galactic disc to the line of 
sight and the degree of order in the magnetic field ( |Stil et al.|2009| . 
Thus, radio continuum surveys that cover both total and polarised 
synchrotron emission, using the SKA and other radio telescopes, 
will be able to produce statistical distributions of fractional polari¬ 
sation for the local universe and for different redshifts, that can be 
directly compared to our (or similar) models. This data will provide 
a new window through which to test both galaxy evolution models 
and theories of magnetic field evolution in galaxies. For example, 
all of the galaxy formation models considered here predict that less 
than 50 per cent of the spiral galaxies with stellar mass greater than 
10 10 Mq should host large-scale magnetic fields greater than 1 gG 
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in strength. Similar statements are made for other mass ranges in 
the main text and illustrated in Fig. [5] 

Our results suggest that the strength of the large-scale galac¬ 
tic magnetic field increases nearly exponentially with redshift in 
all mass ranges. This does not, however, mean that the magnetic 
field in an individual galaxy evolves exponentially because indi¬ 
vidual galaxies migrate between the mass ranges and, more impor¬ 
tantly, individual galaxies become dynamo active at different stages 
of their evolution (even within a relatively narrow mass interval) 
because of different individual merger trees. 

Our approach can be improved in several ways. Most impor¬ 
tantly, by allowing for a finite magnitude of the dynamo time scale, 
we shall be able to include detailed evolution of magnetic field al¬ 
lowing us to predict with confidence the appearance of individual 
galaxies in the radio range and to assess the importance of mag¬ 
netic fields effects on star formation in evolving galaxies, and to 
trace back a given magnetic configuration to earlier stages. 
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APPENDIX A: THE SCALE HEIGHT OF A 
PRESSURE-SUPPORTED THIN DISC 


GALFORM assumes that the stars and gas in galaxy discs have an 
exponential radial profile of the stellar surface density, 




_ M */s -r/r s 

2 nrf 


(Al) 


where r s is the scale length and this is assumed to be the same for 
stars and gas. 

Under the assumption of hydrostatic equilibrium in the z- 
direction, one obtains (using cylindrical coordinates) the following 
expression for the gas pressure, P , 


d ( 1 dP\ 

^r 47rGpt ’ (A2) 

where p t is the total gravitating mass density. Defining the surface 
densities 

rz gy 

Z(R,z)=2j p(R,z)dz, Ijjr == 2p(R,z), (A3) 

which is related to the total surface density through 
I(R) = lim ; ^„E(/?,z) and using equation (| A2|), we obtain 
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(A4) 


Since stars and gas are assumed to be distributed similarly, we have 
Eg = kLi with k a constant at a given time, which leads to 

did, 
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dz 


Then, integrating equation l |A4| l, gas pressure follows ; 

P=^GI g (I g + E*). 

On the other hand, 


(A5) 


P= 3 V5P ’ 


(A6) 


where a (> 1) allows for various contributions to the interstellar 
gas pressure: thermal, turbulent and from cosmic rays and magnetic 
fields. 

The scale height is obtained by combining equations {AT}, 



Magnetic fields and galaxy formation 


( |A5[ l. {A6} and p ~ T,/h, i.e., assuming weak vertical density vari¬ 
ations: 


h(r) = 


av 0 r s 
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APPENDIX B: GALACTIC OUTFLOWS DRIVEN BY 
SUPERNOVAE 


We follow equation (2) from Mac Low & Mc Cray: (1988) for the 
radius, of a superbubble produced by an OB association, 

( L 38 r 7 3 \ 1/5 

R sb = 267 pc I ~^ 7 I , 


where L 38 is the mechanical luminosity in the unit of 10 38 ergs _1 , 
?7 is time in the unit of 10 7 yr and no is the number density of the 
ambient gas. The superbubble expands at the speed 

/? s b = 15.7kms _1 L3g 5 n 0 2 ^ 5 . 

|Mac Low & McCray] ( | 1988| l find that the superbubble breaks out 
from the disc when R ~ 2 h. This can be used with the previous 
equations to obtain the supperbubble age at the breakout, 
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where /7 33 4 = A/134pc. Therefore, the velocity of the shock front 
produced by the superbubble at the break out is 


t'sh - R sb\ Rsb=2 h 


= 4 kms 
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and assuming L 38 = 1, one arrives at equation d- 












